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ABSTRACT 

We use data on the image size of the radio afterglow of GRB 030329 (Taylor et al. 2004) to constrain the physi- 
cal parameters of this explosion. Together with the observed broad band spectrum, this data over-constrains the 
physical parameters, thus enabling to test different GRB jet models for consistency. We consider two extreme 
models for the lateral spreading of the jet: model 1 with relativistic expansion in the local rest frame, and 
model 2 with little lateral expansion as long as the jet is highly relativistic. We find that both models are con- 
sistent with the data for a uniform external medium, while for a stellar wind environment model 1 is consistent 
with the data but model 2 is disfavored by the data. Our derivations can be used to place tighter constraints on 
the dynamics and structure of GRB jets in future afterglows, following a denser monitoring campaign for the 
temporal evolution of their image size. 

Subject headings: gamma-rays: bursts — ISM: jets and outflows — radiation mechanisms: nonthermal — 
polarization — relativity — shock waves 
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1. INTRODUCTION 

It has long been recognized that direct imaging of 
Gamma-Ray Bursts (GRBs) can pr ovide important con- 
straints on their physical parameters ("Granot & Loeb| ; 
lOranot, Piran & Sari 1999a b; Panaitescu & Meszaros 11998 ; 
ISariiil99& Waxman 1997). Unfortunately, the characteris- 
tic size of a GRB image is only of order a micro-arcsecond 
about a day after the GRB at the Hubble distance, and so it 
cannot be resolved by existing telescopes. Nevertheless, indi- 
rect constraints on the image size of GRB afterglows were 
derived based on the transition between diffractive and re- 
fractive scintillations (Goodman 1997) in the radio afterglow 
of GRB 970508 (Frail et al. 1997; Waxman. Kulkarni & Fraill 
[D98), and based on microlensing by a star in a foreground 
galaxy (Loeb & Perna 1998) fo r the optical lig htcurve of GRB 
00030 IC (Garnavich, Loeb & Stanek' '2000; 'Gaudi & Loeb' 
II^Ol; Gaudi. Granot. & Loeb 2001; Granot & Loeb 2001; 
lMao&T,oebl2001D. 

Obviously the challenge of imaging a GRB is made easier 
for nearby sources where the late ra dio afterglow extends over 
a wide, possibly resolvable angle ( Cen 199?;' Granot & Loeb 
12003; Paczvriski 2001; Wang & Loeb 2001; Woods & Loeb 
[1999). Recendv. ITavlor et alj (2004) have used a VLBI cam- 
paign to measure, for the first time, the angular size and proper 
motion of the radio afterglow image of the bright, nearby 
(z = 0. 1685) GRB 030329. The diameter of the afterglow im- 
age was observed to be ^ 0.07 mas (0.2 pc) after 25 days and 
0.17 mas (0.5 pc) after 83 days, indicating an average veloc- 
ity of ~ 4.1-5.7 c. This superluminal expansion is consis- 
tent with expectations of th e standard relativistic jet model 
iOien. Nakar & Piranl2004 V The projected proper motion of 
GRB 030329 was measured to be smaller than 0.3 mas for 80 
days following the GRB. 

Here we use the data of ITavlor et"al to con- 

strain the physical parameters of GRB 030329 based on de- 
tailed modelling of the collimated relativistic hydrodynam- 
ics of GRB afterglows. Since the current state-of-the-art 
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modelling of afterglow jets is still flawed with uncertain- 
ties (Cannizzo. Gehrels & Vishniac 2004; Granot et al. 200 jl 
Kumar & Granot 2003; Rhoadsj il999t iSalmonsoni 120031 
Sari. Piran & Halpern 1999), we use this data to critically 
assess some classes of models that were proposed in the 
literature. An important difference between relativistic ra- 
dio jets of G RBs and the better-studied relativi stic radio 
jets of quasars (Begelman. Blandford & Rees 1984) or micro- 
quasars (Mirabel & Rodriguez 1999) is that active quasars of- 
ten inject energy over extended periods of time into the jet 
while GRB sources are impulsive. Although quasar jets re- 
main highly collimated throughout their lifetimes, GRB jets 
decelerate and expand significantly once they become non- 
relativistic, ^ 1 yr after the explosion. The hydrodynamic 
remnant of a GRB eventua lly becomes nearly spherical only 
after -5x10^ yr (lAval & Piran.200"ll) . 

The outline of the paper is as follows. In ^ we discuss 
the expected image size of radio afterglows and its relation to 
the observed flux density below the self absorption frequency. 
In ^ we analyze the expected temporal evolution of the af- 
terglow image size. The expected linear polarization is dis- 
cussed in 13 while the surface brightness profile across the 
image and its effects on the estimated source size are consid- 
ered in ^ Finally, we apply these derivations to the radio 
data of GRB 030329 (^ and infer the physical parameters 
from the measured spectrum (^. We conclude in |8lwith a 
discussion of our primary results and their implications. 

2. THE IMAGE SIZE AND SYNCHROTRON SELF ABSORPTION 

In GRB afterglows, relativistic electrons are accelerated in 
the advancing shock wave to a power law distribution of en- 
ergies, dN /d^e oc %^ for 7(. > 7^. For p > 2, the minimal 
Lorentz factor of the electrons is given by 



1 



^e.(r-i), 



(1) 



where is the fraction of the internal energy behind the shock 
in relativistic electrons, and F is the bulk Lorentz factor of 
the shocked fluid. There is a spectral break at = J^syn(7m), 
the synchrotron frequency of electrons with 7(. = 7,,, . Another 
break in the spectrum occurs at i^c = i^synijc), the synchrotron 
frequency of an electron that cools on the dynamical time. 
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At sufficiently low frequencies, below the self absorption 
frequency Vsa, the optical depth to synchrotron self absorp- 
tion becomes larger than unity, causing an additional break 
in the spectrum. In this spectral range, the emitted intensity is 
given by the Rayleigh- Jeans part of a black body spectrum, 
where the black body temperature is taken as the effective 
temperature Teff of the electrons that are emitting the radia- 
tion at the observed frequency v. In the local rest frame of the 
emitting fluid this may be written as 



(2) 



where primed quantities are measured in the local rest frame 
of the emitting fluid while un-primed quantities are measured 
in the observer frame (the rest frame of the central source). 
When Vsa > v„,, the emission at v„, < v < v.a is dominated 



by electrons for which v ^ !^syn(7e) oc 7^, giving 7eff oc v 



1/2 

w - ^synv lej ^ le^ ;eii ^ " 

and oc oc v^l'^ . For v„, > there is rapid cooling and 
all the electrons cool significantly within a dynamical time 
dSari. P iran & Naravan 1998). When v„, > max(i'c, i^m), then 
as v decreases below Vsa the distance / behind the shock where 
Tj/(/) = I decreases. The electrons in that location, which 
are responsible for most of the observed emission, have had 
less time to cool after passing the shock and therefore have 
a higher T^ff = 7eff(ffJeC^/^B). In this case jeff oc 1// cx i^"^/** 
and /v oc z^"/*^ ("Gra not. Piran & Sarill2000h . At a sufficiently 
small distance behind the shock, smaller than Ic, an elec- 
tron with an initial Lorentz factor 7,„ does not have enough 
time to cool significantly after crossing the shock. There- 
fore, most electrons within a distance of Ic from the shock 
have 7(, ~ 7^, and the effective temperature in this region is 
r^ff ?^ ^„mpC^ /kn. At sufficien tly low frequencies (below Vac, 
see lGranot. Piran & S ari 2000) / becomes smaller than Ic and 
7eff ~ 7m independent of v, and therefore oc v < Vac- 
For slow cooling {v,n < Vc), 7eff ~ 7m and /v oc immedi- 
ately below Vsa- 

The observed specific intensity is given by Ii, = {vjv')^I'^^, 
and v' jv = (1 +z)r(l - /3cos6') (1 +z)/r where z is the 
source redshift and 6 is the angle between the direction to 
the observer and the velocity vector of the emitting ma- 
terial in the observer frame. The observed flux density 
is iv = JdncosOI^ « ni^ where « Tr(R±/DAf = (1 + 

z)MR± /Dpf = ( 1 + zfiriRs^ /Dif and 6^ ^ tan 6* = /£>a < 
1 are the solid angle and angular radius of the source im- 
age, respectively. Here R±^ is the radius of the observed 
image (its apparent size on the plane of the sky) and Da, 
Dp and Di are the angular, proper and luminosity distances 
to the source, respectively. Thus one obtains I^, [r/(l + 
z]]3[2(zy')Vc2];t7;ff « [r/(l+z)]2iy27gffm„ and jKatz & PiranI 



(3) 



In deriving Eq. the specific intensity was assumed 
to be uniform across the observed image. A more accu- 
rate calculation would have to integrate over the contribu- 
tion to the observed emission from different radii R and an- 
gles 9 from the line-of-sig ht for a fixed observed time t (e.g., 
iGranot. Piran & Sarilll999 b). which results in a non-uniform 
Ii, across the image. Therefore, when using Eq. one must 
choose some effective value for which should correspond 
to its average value across the image. Since depends on F, 
this is equivalent to choosing an effective value of F. Since F 




Fig. 1. — Schematic illustration of the equal-arrival time surface (thick 
black line), namely the surface from where the photons emitted by the shock 
front arrive at the same time to the observer (far on the right-hand-side). The 
maximal lateral extent of the observed image, Rx, is located at an angle 
6t, where the shock radius and Lorentz factor are R» and F, = rsi,(i?»), 
respectively. The area of the image on the plane of the sky is S± = ttR-^ . 
The shock Lorentz factor Fsh varies with R and B along the equal-arrival 
time surface. The maximal radius, R;, on the equal-arrival time surface is 
located along the line-of-sight. If, as expected, Fsh decreases with R, then 
T/ = Fs[,(R/) is the minimal shock Lorentz factor on the equal-arrival time 
surface. 



depends on R, one also has to find at which R 01 9 should 
the value of F be evaluated in Eqs. Q and Q. Usually 
Ksa < i^m < Vc in which case 7eff w 7,,, so that ly depends 
on F not only through the Lorentz transformations, but also 



(|3j. Comp aring E q. (|3 with the more accurate expression cal- 



through the value of_7,„, i.e. F enters into both Eqs. Q and 
t aring E 

culated bv lGrano t & Sari (2002) using the Blandford-McKee 
(1976) self similar spherical solution, we find that the two ex- 
pressions are in relatively good agreement'* if F is evaluated 
just behind the shock at the location where R±^ is located. This 
should be a good approximation before the jet break time in 
the light curve. 
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1/(3-*) 
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'0.66(1 +z){E5i/noy^\9o/Q.lf days (k = 0) 

0.34(1 +z)iE5i/A,)i9o/Q.lf days (k = 2) 

At f > tj, however, it is less clear how well this approximation 
holds, and it might be necessary to evaluate F at a different 
location. In particular, as we shall see below, this approxima- 
tion does not work well for model 2 with k = where F need 
to be evaluated near the head of the jet, rather than at the side 
of the jet where R± is located. 

The image size is given hy R± = max(Rsm9) along the 
equal arrival time surface (see Figure The equal arrival 
time surface is the surface from where photons that are emit- 
ted at the shock front arrive to the the observer simultane- 
ously. Since the emission originates only from behind the 
shock front, the projection of the equal arrival time surface 
onto the plane of the sky (i.e. the plane perpendicular to the 
line-of-sight) determines the boundaries of the observed im- 
age, and its apparent size (see Figure[3. For a spherical shock 
front with any Ritj^b), R± = inax(Rsm9) is located at an an- 
gle 6** which satisfies cos 6** = /3* (see Appendix IaI. where 
and F* = (1 - are the velocity (in units of c) and 

the Lorentz factor of the shock front^ at 6^ . This implies that 
R±{t) = RJt) /F*(f ) where RJt) = R(t, cos 9 = [3^) is the radius 

^ We find that the ratio of the numerical coefficient in Eq. |3) to that 
in Granot & Sari i2002) is in this case L09(3p- l)/(3p + 2) ior k = 2 and 
L71(3p-l)/(3p + 2) for/t = 0. 

^ Note that we use fti, or Fjt, for the location of the emitting fluid, which 
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Fig. 2. — Evolution of the jet half-opening angle 0, as a function of radius 
R, for various illustrative cases. The solid line shows the evolution derived 
from 2D hydrodynamical simulations I Granot et al. 2001 ). The different lines 
give the maximal polar angle d of the shock front (which is obtained at a rel- 
atively small radius where a minor fraction of the emission is produced), and 
the average values of d within the jet when averaged over the circumburst 
gas density and over the total emissivity. Most of the emission comes from 
within the original jet opening angle, = 0-2. Also shown is the evolu- 
tion oi 9j{R) predicted by simple semi-analytic models. Thi'ee illustrative 
cases are depicted where the lateral expansion speed is assumed to be /3j = 0, 
3-'/^ and 1 in the local rejt^frame |^en. Nakai" & Piran 2004; Rhoads 1991 
ISari. Piran & H alpern 1999). Since the onset of lateral expansion in the sim- 
ple models takes place at a somewhat larger radius (Rj), a higher value of 
gas density is adopted for these models in order to show more easily their 
different qualitative behaviors. 



of the shock at = arccos/3». Therefore = FshC/?*) and 
P* = PshiR*)- Although the shock front is probably not sim- 
ply a section of a sphere ( Granot et al. 2001 ), we consider this 
as a reasonable approximation for our purpose. The expres- 
sion for 6'* in the more general case of an axially symmetric 
shock is given in AppendixlAl 

The apparent speed, /3ap = [(l+z)/c]{dR±/dt), has a simple 
form for a point source moving at an angle 6 from our line-of- 
sight, /3ap = /3shsin6'/(l -/3shcos0). Substituting cos9 = fish in 

this expression gives /3ap = r,/3* = \/Tl - 1 or = ^Jl + 0l^ 

and /3h. = Ap/ \J^+Pl^- In Appendix 151 we show that this 

result holds for any spherically symmetric shock front, and 
we also generalize it to an axially symmetric shock. Fi- 
nally, the Lorentz factor F of the shocked fluid just behind 
the shock at 0* is related to the Lorentz factor of the shock 
itself, F., by F; = (F+ l)r^ (F- 1)+ llVr^(2-7)(F- l) + 2] 
jBlandford & McKeel 19761) where 7 is the adiabatic index of 
the shocked fluid. For F* > 1, 7 = 4/3 and F = F*/\/2. 

3. THE TEMPORAL EVOLUTION OF THE IMAGE SIZE 

For simplicity, we consider a uniform GRB jet with sharp 
edges and a half-opening angle Oj, with an initial value of 6*0. 
The evolution of the angular size of the image and its angular 
displacement from the central source on the plane of the sky, 
for viewing angles 6'obs > 6*0 from the jet axis, was outlined 
in Granot & Loeb (2003). Here we expand this discussion to 
include viewing angles within the initial jet opening angle, 
^^obs < Sq, for which there is a detectable prompt gamma-ray 
emission (similarly to GRB 030329 which is considered in 
the next section). For 6'obs < 6*0, Ri^ is the observed size of the 

is always just behind the shock. On the other hand, we use /3 or F (which are 
slightly smaller than /3sh or V^^ , respectively) for the Lorentz transformations 
of the emitted radiation, since these are the bulk velocity and Lorentz factor 
of the emitting fluid. 



image, while for 0obs > ^'o it represents the displacement with 
respect to the central source on the plane of the sky. 

In this section we concentrate on a viewing angle along 
the jet axis, 6'obs = 0, and in the next section we briefly out- 
line the expected differences for < 6'obs < ^'o- For ^obs = 0, 
the observed image is symmetric around the line-of-sight (to 
the extent that the jet is axisymmetric). At f < tj the edge 
of the jet is not visible and the observed image is the same 
as for a spherical flow: cx (£i,o/A)i/2(4-*)f(5-/:)/2(4-/:) ^ 
(£/A)i/2(3-<:)f-i/2(4-t)j(5-i:)/2(4-i:) f^j. gxtemal density profile 

pg^t = Ar"*, i.e. a = {5-k)/2{A-k) where a = d\\\R^/d\nt. 
Here E is the true kinetic energy of the jet, and Ei^o = f^^E 
is the isotropic equivalent energy where = 1 - cos 6q w 
6'q/2 is the beaming factor. At f < tj the flow is de- 
scribed by the Blandford-McKee (1976) self similar solu- 
tion, wh ich provides an accurate expression for the image size 
JGranot. Piran & Sari 1999a: Granot & Sari 2002,) . 

^^~[ 7r(5-/fe)5-*(l+z)5-':A 

r 3.91 X 10i«(£„o.52/«o)'/''[fdays/(l+z)]'/''cm (k = 0) 

= { (5) 
[2.39X 10i'>(£iso,52/A,)'/4[f,,y,/(l+z)]3/4cm (k = 2) 

At f > the jet gradually approaches the Sedov-Taylor self 
similar solution, asymptotically reaching R±^ cx (£'f^/A)'/'^^~*', 

1. e. a = 2/(5-k). At tj < t < Jnr there is a large uncertainty 
in the hydrodynamical evolution of the jet, and in particular 
its rate of sideways expansion. We therefore consider two ex- 
treme assumptions which should roughly bracket the different 
possible evolutions of R \ (t): (1) re lativistic lateral expansion 
in the comoving frame (' Rhoadslll999t ISari. Piran & Halpernl 
1999), for which dj w max(6'o,7"^) so that at tj <t < ?nr we 
have 7 « 61"' w 6'q' exp(-R/Rj), and (2) little or no lateral ex- 
pansion, 6j w 6() for t < fNR, in which case appreciable lateral 
expansion occurs only when the jet becomes sub-relativistic 
and gradually approaches spherical symmetry. We shall refer 
to these models as model 1 and model 2, respectively. Model 
2 is also motivated by the results of numerical simulations (see 
Figure 13 which show only modest lateral expansion as long 
as the iet is relatiy istic ( Cannizzo. Gehrels & Vishniac 2005 
iGranot et all 1200 U iKumai- & Granot 2003). These numeri- 
cal results are also supported by a simple analytic argument 
that relies o n the shock jum p conditions for oblique relativis- 
tic shocks ( Kumar & GranoL20 03). 

Figure |3] schematically shows the evolution of R±{t) 
for these two extreme models, both when viewed on-axis 
(6'obs < ^o) as required for seeing the prompt gamma- 
ray emission, and for 0obs ~ 90° as will typically be 
the case for GRB jets found in nearby SNe Type 
Ib/c fc ranot & Loe b l2003t iG ranot & Ramirez-Ruiz 2004 
lPaczvhskiiii2001i iRamirez-Ruiz & Madau..200 4). For 6'obs < 
610 at < f < fNR we have R± oc (£/A)>/2(3-/:)f 1/2 f^j. j^j^^j^j ^ 
andR± cx (£i,sof/A)V(4-*) ^ (£/A)i/(3-i)(f/f .)i/(4-« for model 

2. Therefore, with k = 2 we have a = 1 /2 for both models, 
despite their very different jet dynamics. For A: = we have 
fl = 1/2 for model 1 and a = 1/4 for model 2. 

In model 1, the jet is already relatively close to being spher- 
ical (i.e. dj 1) at fNR = fNR(£'), where R^r{E) = cInrIe) = 
[(3 - A:)£'/47rAc^] and its radius is similar to that of 
the Sedov-Taylor solution , RsT(E,t) = .^(EfVA)'/^^"**, cor- 
responding to the same time f, where ^ = ^(^,7) ^ 1. There- 
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Fig. 3. — Schematic plot of the evolution of the observed afterglow image 
size R± of a uniform GRB jet with sharp edges. The jet is either viewed from 
within the initial jet opening angle, Sqi,, < 8q (upper panel), or from d^bs ~ 
90° (lower panel). The solid line is for model 1 (relativistic lateral expansion 
in the local rest frame) and the dashed line is for model 2 (little or no lateral 
expansion before ^nr). The dotted line in the upper panel represents jets (in 
model 2) with a smaller dg and the same true energy E, which converge to the 
same self similar dynamics and therefore the same R±(t) after the jet break 
time tj. Also shown are the ratios of various values of and t. 



fore, we expect it to approach spherical symmetry on a few 
dynamical times, i.e. when the radius increases by a factor of 
b ^ Si few, corresponding to a factor of ^ b'-^~'^^l^ in time, and 
the transition between the asymptotic power laws in R±_{t) is 
expected to be smooth and monotonic. 

In model 2, however, the jet becomes sub-relativistic 
only at /?NR(£iso) = ct^RiEiso), which is a factor of ~ 
(£i.o/£)'/(3-*) = fl'"'-"' ^ e-^"'-'' larger than R^^{E) = 

ctj^R(E) and a factor of ^ fi,^^''^ ^ ^o^^*^ ''^ larger than 
/?ST[£^,fNR(£iso)]- It also keeps its original opening angle, 
9j ^ 9o until fNR(£'iso), and hence at this time the jet is still 
very far from being spherical. Thus, once the jet becomes 
sub-relativistic, we expect it to expand sideways significantly, 
and become roughly spherical only when it has increased its 
radius by a factor of b ^ a few. This should occur, however, 
roughly at a time r^ph when R5T(E,t^j,h) = MNR(£'i,so), i-e. 



^sph 



V20,7'/7(5-*)/2 . (6) 



This is a factor of ^ '"'^ w 14(0()/O-1) ' larger than the ex- 
pected transition time in model 1, and for b ^ 2-3 gives 
a factor of - (80 - 220)(6'o/0.iri for ;t = and - (40- 
70)(6'o/0.1)"' forA: = 2. During this transition time, /?j_(6'obs < 
6*0) grows by a factor of - f'^^^b - d^^^b while R±{9obs ~ 90°) 
grows by a factor of ^ b. This implies that during the transi- 



lnMl/2)ln/i, 
[(5-i)/2]lnMl/2)lnA 

Inb 

[(5-i)/2]lni-(l/2)ln/6 



(6'obs < ^o) 



(^obs = 90°) 



(7) 



and < (fl) < 2/(5-k) for ^obs < 6*0 while 2/(5 -/t) < (a) < 
1 for 6'obs = 90°, where (a) « 2/(5 -/t) in the Hmit b > 6*0' 
(which is not very realistic). The other limiting value of (a) w 
for 6'obs < ^'o and (a) w 1 for 0obs = 90° is approached in the 
limit b<^9^^^ . Typical parameter values (b ^ 2-3, 9o ^ 0.05- 
0.2) are somewhat closer to the latter limit. For example, for 
/t = 0, = 2.5 and 9o = 0. 1 we have (a) w 0.722 for 6'obs < 6*0 
and (a) w 0. 185 for 6'obs = 90°. This demonstrates that for on- 
axis observers there should be a sharp rise in R±, while for 
observers at 6'obs ~ 90° there should be a very moderate rise 
in R± during the transition phase from the asymptotic tj <C 
t <C fNR and t ^ fNR regimes. Furthermore, as is illustrated in 
Figure|3l this transition would not be monotonic in model 2. 
This is because during the transition a passes through values 
larger (smaller) than both of its asymptotic values for 6'obs < 9o 
(0obs « 90°). 

For comparison, and in order to perform a quantitative com- 
parison with the data, we consider a simple semi-analytic 
model where the shock front at any given lab frame time 
occupies a section of a sphere within 9 < 9j, and R± is lo- 
cated at 9± = mm{9^,,9j). The observer time assigned to a 
given 6'_L(fiab) is t = fiab - [^(fiab)/c] cos6'^(fiab). We follow 
Oren. Nakar & Piran (2004) with minor differences: (/) we 
choose the normalization of R± at f <C tj so that it will co- 
incide with the value given by the Blanford-McKee solution 
(i.e. Eq. and (ii) the lateral spreading verlocity in the co- 
moving frame, Pj, for model 2 smoothly varies from f3j <C 1 
at f <C fNR to the sound speed, jSj « q/c, at f > Jnr. The latter 
is achieved by taking Pj to be the sound speed suppressed by 
some power of F. 

Figurel^shows the resulting R±{t) for ISM (k = 0) and stel- 
lar wind {k = 2) environments, and different recipes for f3j. For 
a given (3j recipe, R±{t) depends on E/A and 6'o. The values 
of these parameters that were used in Figure |3 are indicated 
in the figure. For k = 2 the spread in R±{t) for the different 
(3j recipes is smaller than for k = Q. This is understandable 
since the asymptotic values of a are the same for models 1 
and 2. There is still a non-negligible spread, however, as the 
asymptotic value of a = 1 /2 at f ^ <C f <C fNR is not reached.^ At 
t ^ fNR all recipes for /3j approach the same value of R±{t), 
except for f3j = for which ^_L(f) is smaller by a factor of 
sin 00- For (3j = and k = there is a pronounced flattening 
in R±{t) at ^ 1 .2 day, which is a factor of ^ 7 larger than the 
value of tj = 0.165 days that is implied by Eq. (|4j. We must 
stress that this simple model becomes unrealistic around fNR. 

The apparent velocity of a point source is /?ap = /?sin6'/(l - 
(3cos9). For 6'obs > 9o, as long as 9j < 9obs and f < fNR we 
have /3ap w 2r2j,6l/[l + (rsh6l)2] « 2/61. For 6lobs = 7r/2 we have 
/3ap = Psh which is close to 1 at f < fNR. For Fsh ^ 1 and 
9 > r^jj we have /3ap ~ sin6'/(l - cos 6'), so that /?ap > 1 for 
^obs < 7r/2 and /?ap < 1 for 6'obs > tt/2 (i.e. fo r the counter jet 
assuming a double sided jet; see Figure 2 of iGranot & Loehl 
2003). For 6'obs < ^'o we have /3ap = r*/3,^ at f < fNR. At 
f < tj we get 6'j^ = 6'* < 6'() and the shock front is roughly spher- 
ical with an approximately uniform Lorentz factor within 

* This is since it takes a long time to approach this limit for k = 2, which is 
longer than the dynamical range between tj and Jnr. 
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Fig. 4. — Evolution of the source size (or more precicely, its diameter 2R±) 
as a function of time, for a uniform density environment (k = 0, upper panel) 
and for a stellar wind (k = 2, lower panel). Different recipes are considered 
for the lateral spreading velocity in the comoving frame, f3j. See text for 
more details. 



9<9^,,so that T^, w T^h- At tj<t < we have 9^ « 9j w r;:jj 
and /3ap w w Fsh for model 1, suggesting that using r(6'*) 
for calculating the emission (i.e. in Eqs. ^ |3] and |8|i is a 
reasonable approximation. For model 2, 6'* « 9j w 9q and 
/?ap « F,, w 29oTl^^, so that^ F^/F^h w 26ioFsh < 1 suggesting 
that F(6'*) underestimates the effective value of the emissivity- 
weighted F , which enters the expressions for the observed 
emission. This results from the fact that in model 2 most of 
the emission at tj <t < Jnr originates fr om 9 <9n where F is 
higher than at 9^ > 9o (see Figurel2land lGranot et alJ20Qll) . 

4. LINEAR POLARIZATION 

For < ^obs < ^'o the image would not be symmetric around 
the line-of-sight, but its typical angular size would be simi- 
lar to that of 6'ob.s = 0. If there is significant lateral spreading 
at f > tj, then this should cause the image to become more 
symmetric around our line-of-sight with time. This, by itself, 
might be a possible diagnostic for the degree of lateral spread- 
ing. The degree of asymmetry in the observed image should 
also be reflected in the degree of linear polarization, and its 
temporal evolution. While the image might be resolved only 
for a very small number of sufficiently nearby GRBs, the lin- 
ear polarization might be measured for a larger fraction of 
GRBs. 

Contrary to naive expectations, for very slow lateral ex- 

' Here Fsj, represents the uniform shock Lorentz factor in the simple semi- 
analytic model described at the end of ^ where the shock at any given fiab 
occupies a section of a sphere and abruptly ends at 9j, and at tj < t < fNR 
R± is located at 9j. On the other hand, F, = FshC^*) is the Lorentz factor at 
the angle 6» where is located for a smooth and continuous (and therefore 
more realistic) shock front, for which changes with 6 at a given fi^b ■ 



pansion (J3j ^1) the polarization decays faster after its peak 
at f ~ tj compared to lateral expansion at the local sound 
speed, l3j = cjc « 3"'/^, in the comoving frame ("Rossi et alj 
^(34). A very fast lateral expansion in the local frame close 
to the speed of light {(3j « 1), leads to 9j k, max( 6*0,7"') 
and to three peaks in the polarization light curve, where the 
polarization position angle changes by 90° as the d egree 
of polarization passes through zero between the peaks JSari 
[1999). When there is a slower lateral expansion or no lat- 
eral expansion at all ( Ghisell ini & Lazz ati 1999), there are 
only two peaks in the polarization lightcurve where again 
the polarization position angle changes by 90° as the degree 
of polarization passes through zero between the peaks. The 
peak polarization is higher for f3j w (~ 15%- 16%) com- 
pai-ed to I3j = 3"'/^ (- 9%) jRossi et alJl2004h . The maximal 
observed degree of polarization is, however, usually < 3% 
suggesting that the magnetic field configuration behind the 
shock is more isotropic than a random field fully within the 
plane of the shock (Granot & Konigl 2003) which is expected 
if the magnetic field is produced by the Weibel instability 
jMedve dev & Loeb 1999). This changes the overall normal- 
ization of the polarization light curve, and hardly affects its 
shape (Granot & Konigl 2003). Since the overall normaliza- 
tion is the most pronounced difference between slow and fast 
lateral expansion, and it is very similar to the effect of the de- 
gree of anisotropy of the magnetic field behind the shock, it 
would be very hard to constrain the degree of lateral expan- 
sion from the polarization light curves. There are also other 
possible complications, such as a small ordered magnetic field 
component (Granot & Konigl 2003) which can induce polar- 
iz ation that is not related to the jet structure. 

iTavlor et' al. (2004) put a 3 ct upper limit of 1% on the linear 
polarization in the radio (i^ = 8.4 GHz) at f = 7.71 days. They 
attribute the low polarization to synchrotron self absorption. 
Indeed, i^sa is above 8.4 GHz at this time, but only by a factor 
of ^ 2. One might expect a suppression of the polarization in 
the self absorbed region of the synchrotron spectrum since it 
should follow the Rayleigh-Jeans part of a black body spec- 
trum, and depend only on the electron distribution (i.e. the 
"effective t emperature") and not on th e details of the mag- 
netic field ( Granot. Piran & Saril999bh . The optical depth to 
self absorption does, however, depend on the details of the 
magnetic field, and may thereby vary with the direction of 
polarization. Therefore, there might still be polarization at 

^ t^sa which will go to zero in the limit ly <C Kso- An or- 
dered magnetic field in the shocked fluid through which the 
emitted synchrotron radiation propagates on its way to the ob- 
server, might induce some polarization in the observed radi- 
ation (Sagiv, Waxman & Loeb 2004). These effects are sup- 
pressed roughly by a factor of the square root of the ratio be- 
tween the magnetic field coherence length and the width of 
the emitting region (which is of the order of the typical path 
length of an emitted photon through the shocked plasma be- 
fore it escapes the system). 

5. THE SURFACE BRIGHTNESS PROFILE 

ITavlor et all (12001 use a circular Gaussian profile for their 
quoted values, and also tried a uniform disk and thin ring. 
They find that a Gaussian with an angular diameter size of 
1 mas is equivalent to a uniform disk with an angular diam- 
eter of 1 .6 mas and a thin ring with an angular diameter of 
1 . 1 mas. At f < tj the jet dynamics are close to that of a spheri- 
cal flow, since the center of the jet is not in causal contact with 
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its edge, and the dynamics can be described by the Blandford- 
Mckee (1976) spherical self similar solution (within the jet, 
at 9 < 9o). The surface brightness in this case has been in- 
vestigated at length in severa l works ( Granot & Loeb 200 u 
lyranot. Piran & Sari 1999a b; Panaitescu & Meszaros 199g 
ISari 1998 ; Waxman 1997). The surface brightness profile 
of the image, normalized by its average value across the 
image, is the same within each power law segment of the 
spectrum, but changes between different power law segments 
dGranot & Loeb 2001). The afterglow image is limb bright- 
ened, resembling a ring, in the optically thin part of the spec- 
trum and more uniform, resembling a disk, at the self ab- 
sorbed part of the spectrum. This can affect the angular size of 
the im age that is inferred from the observations ( Tavlor et ^ 
iooi, where the angular diameter for a uniform disk (thin 
ring) is a factor of 1.6 (1.1) larger than the values quoted by 
l^vlor et al. (2004) for a circular Gaussian surface brightness 
profile. This effect would be more important at < lyga where 
the afterglow image resembles a uniform disk rather than a 
thin ring. 

One should keep in mind that the image size of GRB 
030329 was inferred well after the jet break time, f ^ tj, and 
relatively close to the non-relativistic transition time, f ^ Jnr. 
However, at tj <t < Jnr the jet dynamics is poorly known, and 
this uncertainty must necessarily be reflected in any calcula- 
tion of the afterglow image at this stage, which could only be 
as good as the assumed dynamical model fo r the jet. The after- 
glow image at this stage was calculated bv lloka & Nakamural 
12001) assuming lateral expansion at the local sound speed 
(iRhoads 1999), similar to our model 1. They find that at t < tj 
the surface brightness diverges at the outer edge of the image, 
which is an artifact of their assumption of emission from a two 
dimensional surface ( Granot & Loeb,2001: Sari. 1998.) identi- 
fied with the shock front. Calculating the contribution from 
all the volume of the emitting fluid behind the shock makes 
this divergence go away, except for certain power law seg- 
ments of the spectrum where the emissi on indeed arises from 
a very thin layer just behind the shock ( Granot & Loebl200ll) . 
At f > tj loka & Nakamura ( 2001 ) obtain a relatively uniform 
surface brightness profile. However, this is due to the unphys- 
ical assumption that the shock front at any given lab frame 
time is part of a sphere within some finite angle 6j from the 
jet symmetry axis where the jet ends abruptly. The edge of 
the image in this case corresponds to this un-physical point 
where the jet ends abruptly (see Figure |5}. More physically, 
as is shown by numerical simulations ( Granot et al. 200 1 ), the 
shock front is not a section of a sphere and is instead round 
without any sharp edges. Similarly to the spherical-like evo- 
lution at f < tj, the edge of the image would in this case corre- 
spond to Ri^ = max(/? sin 6*), and thus the image is expected to 
be limb brightened for the same qualitative reasons that apply 
at f < tj, even though there would be some quantitative differ- 
ences. A proper calculation of the afterglow image at f > tj 
requires fuU numerical simulations of the jet dynamics. 

6. APPLICATION TO GRB 030329 

We now apply the expressions derived in the previous sec- 
tion to GRB 030329 which occurred at a redshift of z = 
0.1685. We use the image angular diameter size of 9^ w 

70 juas for^ w 587 Mpc that was inferred ait = 24.5 days 
dTavlor et al. 2004) . which corresponds to R±^ w 0.1 pc. This 

^ Throughout this paper we assume Qm = 0.27, JIa = 0.73 and Ho = 

71 kms"' Mpc"'. 
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Fig. 5. — Evolution of the source size as a function of time for t > tj. In 
model 1, the lateral expansion in the local frame is relativistic while in model 
2 there is no lateral expansion at f > tj . The shock front at any given time t^ja 
is assumed to be part of a sphere which abruptly ends at a finite angle 9j from 
the jet axis. The gray dotted lines represent the equal aiTival time surfaces at 
three different observed times. Since the jet dynamics, R(t\ih), are different 
for models 1 and 2, the equal arrival time surfaces should be different (but 
in this sketch, for simplicity, we depicted them as being equal). At ? > tj 
(where 6j(Rt.) < 6* ), the edge of the image which determines the image size 
is located at the edge of the jet, i.e. at an angle Bj instead of 



implies an average apparent speed of (/3ap) = (1 +z)R±/ct sa 
5.66. The instantaneous apparent speed is given by /3ap = 
[(l+z)/c]dRjL/dt = a{l3^p) where a = c/lnTJ^/c/lnf. ForGRB 
030329, if we also take into account the inferred source size 
of 9s ~ 170 /ias or R± w 0.25 pc at f = 83.3 days and the 2 <t 
upper limit of 0, < 100 /ias or R± < 0.14 pc at f = 51.3 days 
dTavlor et al.l2004 . we have** a = 0.71!°^ (1 a). This value is 
between t = 24.5 days and 83.3 days, assuming that R±(t) fol- 
lowed a perfect power law behavior cx t" with a = const during 
this time. This is a reasonable approximation for model 1 or 
model 2 with k = 2 for which a = 1 /2 at < f < Tnr and there- 
fore these models are consistent with the observed temporal 
evolution of the image size. For model 1 with ^ = (see ^ 
a = 1 /4 at ^ f <C fNR but its value is expected to increase 
significantly near Jnr which we find to be at ^ 200 days for 
this model (see Table 1). Therefore, it can still account for the 
observed image size at f = 24.5 days and 83.3 days together 
with the upper limits at 51.3 days. At f = 24.5 days, however, 
we still expect the value of a in model 2 with A: = to be rela- 
tively close to its asymptotic value of a = 1/4. 

Figure|6lshows crude fits between the simple semi-analytic 
realization of models 1 and 2 (that is described at th e end 
of ^ and the observed image size dTavlor et aljl200l . For 
model 2 we have used the recipe 9j = r~'(cj/c) for the lat- 
eral expansion. We have treated the value of E/A as a free 
parameter whose value was varied in order to get a good fit, 
while the value of 9o was determined according to the ob- 
served jet break time tj w 0.5 days using Eq. In the latter 
procedure we take into account an increase in energy by a fac- 

^ A pplying the Bayesian inference formalism developed bv lReichart et al. I 
120011) . we determine values and uncertainties for the model parameter a. 
Bayesian inference formalism deals only with measurements with Gaussian 
error distributions, not with lower or upper limits. However, this formalism 
can be straightforwardly generalized to deal with limits as well, using two 
facts: (1) a limit can be given by the convolution of a Gaussian distribution 
and a Heaviside function; and (2) convolution is associative. 
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Fig. 6. — A tentative fit of a simple semi-analytic realization of models 1 
and 2 to the observed image size (of diameter 2R±). The physical parameters 
and external density profile for each model are indicated. 



tor of / ~ 10 due to refreshed shocks CGranot. Nakar & Piran' 
l2003h between f y and the times when the image size was mea- 
sured. For simplicity, we do not include the effect of the en- 
ergy injection on the early image size. The image size that is 
calculated in this way to not valid before the end of the en- 
ergy injection episode (after several days), but it should be 
reasonably accurate at f > 25 days when its value had been 
measured. The values of E /A and r(25 d) from these fits are 
indicated in Figure|6land in Table 1. 

Fora = (0.25, 0.5, 0.75) we obtain /3apW (1.4,2.8, 4.2), 
(1.7,3.0,4.4) and T(0^) « (1.5,2.4,3.4). The values of r(6l,) 
are similar to the value of F that were obtained from the fit 
to the observed image size for model 1 and model 2 with 
k = 2, but it is smaller for model 1 with ^ = 0, as expected (see 
discussion at end of 

Using the radio data from lBerger et al. (2003). we find that 
iv ~ 10 mJy at f w 25 days and ly = 4.86 GHz which according 
to the spectrum at this time is below Vsg. Berger et al. (2003) 
also estimated the break frequencies at f « 10 days to be i^^a ~ 
19 GHz and i^,,, « 43 GHz, which is consistent with i/ < Vsa < 
Vm at f = 24.5 days. A value of p = 2.25 was inferred for GRB 
030329 (IWillingale et al.ll2004i) . For the power law segment 
of the spec trum where oc i^^ (labeled 'B' in Figure 1 of 
iGranot & S ariii2002i) we have 7eff w 7,,, for which Eqs. ([Q 
and (|3} imply 



1 

2^ 



p-l\ (l+z) 



P-2J r(r-i) 



Da 



(8) 



Using the above values for the flux, R±, r(0*) and p 
for GRB 030329, Eq. (g} gives « (0.10,0.023,0.0099) 
for a = (0.25,0.5,0.75). These values of are some- 
what on the low side compared to the values inferred from 
broad ban d afterglow modelling of other afterglows (e.g., 
iPanaitescu & Kumar 2001b). In Table 1 we show in addition 
to these values of e^, also the values that are obtained when 
evaluating F from the fit to the image size that is shown in 
Figure |6l The largest difference between these two estimates 
of is for model 1 with k = 0, for which evaluating F from 
the fit to the observed source size probably provides a more 
accurate estimate. 

Since Eq. (|8j relies on a small number of assumptions, 
it is rather robust. However, the value of e,, in equation (|8} 
is very sensitive to the value of R±. This is because oc 
l//?ir(r-l)andforF> 1 we have F w r»/%/2 w /3ap/V2 = 



fl(/3ap) /v 2 cx R± so that cx R~^. For example, 6^ = 45 /ias 
(Rj_ = 0.064 pc) at t = 24.5 days, which is stifl within the 
measurement errors, would imply = (0.61,0.14,0.060) for 
a = (0.25,0.5,0.75). The latter values , especially for a » 0.5 , 
are consistent with the value found bv'Willin gale et all ll2004h 
from a broad band fit to the afterglow data: = 0.24 and 
0.18 < ec < 0.31 at the 90% confidence level, and with flie 
value of w 0. 16 found bv lBerger et aP (l2003h . 

7. INFERRING THE PHYSICAL PARAMETERS FROM A SNAPSHOT 
SPECTRUM AT f; < r < /„ 

For model 1, we obtain expressions for the peak flux and 
break frequencies at f ,■ < f < ^nr b y using the expressions for 



t < ti from Granot & Saril 



in order to estimate their 



values at tj, and then using the their temporal scalings at 
t j <t < fNR from Rhoads ( 199^ and|Sari, Piran & Haber^ 
il999). In Appendix Icl we provide expressions for the peak 
flux and break frequencies as a function of the physical pa- 
rameters and solve them for the physical parameter as a func- 
tion of the peak flux and break frequencies for both models 1 
and 2. The results fo r GR B 030329 are given below. 

For GRB 030329. iBerger et al. (2003) infer v,,, ~ 19 GHz, 
Vm ~ 43 GHz and F^, max ~ 96 mJy at t ^10 days, as wel l 
as p = 2.2. Using Eqs. 4.13-4.16 of LSari & EsinI ll2()0l . 



IBer ger et all ll2()0l find £1,0,52 « 0.56i/, 



1/4 
c,13- 



no : 



1. 



< 3/4 



O.lOz^^^j"*, £(, ~ 0.16i'^''[3, and using a value of Oj « 0.3 at 
this time they inferred £51 w 0.25. For the same values of 
the spectral parameters and using our model 1 we obtain 



1/4 



£150,52 = 0.16^-. 

1/4 



1/4 
c,13' 



£51 



: Q36i/Jy-^, no ■■ 



-5/4 
e,13 ' 



= Q.lAvJ, , for ^ = and E\. 



o,52 = 0.10l/,'(t'E51=0.43l/f/^ 



A, = 1.41/^13,6^ = 0.0341/, 



-5/4 
c,13 ' 



= 0.361/^13 for /fe = 



13' 

2. The im- 



plied values ofE/A are shown in Table 1 . The diff erences be- 
tween our values and those of'Berg er et alJ ll2003l) arise from 
differences by factors of order unity between the coefficients 
in the expressions for the peak flux and break frequencies. 
This typically results in differences by factors of order unity 
in the inferred values of the physical parameters. The differ- 
ence in the external density n is relatively large since it con- 
tains high powers of Vsa and v„i (Granot, Piran & Sari I999J 
making it more sensitive to the exact theoretical expressions 
and observational values of these frequencies. 

For model 1 and ^ = we obtain £51 /no = 0.0.024z/^ 



-3/8 

c,13 



compared to £51 /no = 0.1 4i/^ 13 of iBerger etai] ( l200l and 
£51 /no ^ 0.8 that we obtain from the fit to the observed im- 
age size (Figure |6}. Because of the large uncertainty in the 
value of n that is determined from the snapshot spectrum, and 
the large uncertainty in the value of E/n from the fit to the 
image size, these values are consistent with each other within 
their reasonable errors (see Table 1). For model 1 and k = 2 

—1/8 

we obtain £5 1 /A * = . 3 1 i/^, 1 3 compared to £5 1 /A * « . 8 from 
the fit to the observed image size. Here the difference between 
the two values is smaller, but the uncertainty on the two val- 
ues is also smaller (see Table 1). Altogether, the two values 
are still consistent within their estimated errors. 

For our model 2 involving a jet with no significant lat- 
eral spreading, the peak flux is suppressed by a factor of 
(f/f;)-(^-*>/*'*-*' where tj w 0.5 days and t/tj « 20, i.e. a factor 
of « 0.1 1 for k = Q and « 0.22 for k = 2. This implies (see 



appendixO £'iso,52 = 4.7i/^. 13, £51 = 0.21;/^ 13, no = 0.53j/, 
ee = O.Slv'^i*, te = 0.078:/''!3 for ;t = and £1,0.52 = OMiy, 



,3/4 
-.13' 
1/4 

c,13' 
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£51 = 0.291^1^1, A, = lAiy'{l, = 0.071i/^ j^^ e, = O.lVt^'g 

for k = 2. For model 2 with = we get Esi/no = OAOv'^^-^ 
compared to £51 /no ~ 5 from the fit to the observed image 
size. These two values are consistent within the large uncer- 
tainties on both values (see Table 1). 

— 1/8 

For model 2 with k = 2 we obtain £51 /A^ — 0. 10i>'^ j-^ com- 
pared to £51 /A, « 1 .2 from the fit to the observed image size. 
In this case, however, the errors on these two values are rel- 
atively small (see Table 1). This is because: (i) the image 
size is linear in E /A which corresponds to a relatively strong 
dependence, and therefore the observed image size can con- 
strain the value of E /A relatively well, and ( ii) the expression 
for E I A from the spectrum contains relatively small powers 
of the break frequencies and peak flux and thus has a corre- 
spondingly small uncertainty. Therefore, the two values of 
£51 /A* are farther apart than is expected from the uncertainty 
on these values. Thus, one might say that the data disfavors 
model 2 with k = 2. It is hard, however, to rule out this model 
altogether, because of the uncertainty is the exact expressions 
for the break frequencies and peak flux at tj < t < Jnr- 

8. DISCUSSION 

We have analyzed the data on the time-dependent image 
size of the radio afterglow of GRB 030329 (Taylor et al. 
2004) and constrained the physical parameters of this explo- 
sion. The image size was measured after the jet break time tj 
in the afterglow lightcurve, where existing theoretical models 
still have a high level of uncertainty regarding the jet dynam- 
ics. This motivated us to consider two extreme models for 
the lateral expansion of the jet: model 1, where there is rel- 
ativistic lateral expansion in the local rest frame of the jet at 
tj < t < fNR, and model 2, with no significant lateral expan- 
sion until the transition time to a non-relativistic expansion 
fNR- We have tested the predictions of these models against 
the observations, for both a uniform (pext = Ar~*, with k = Q) 
and a stellar wind (k = 2) external density profile. 

The observational constraints included comparisons be- 
tween: ( i) the value of the post-shock energy fraction in rel- 
ativistic electrons that is inferred from the source size and 
flux below the self absorption frequency and its value from the 
'snapshot' spectrum at f « 10 days; ( ii) the value of E /A that 
is inferred from the source size and its value from the 'snap- 
shot' spectrum at f w 10 days; and (Hi) the observed temporal 
evolution of the source size and the theoretical predictions. 



We have found that most models pass all these tests. The only 
exception is model 2 with k = 2, involving a relativistic jet 
with little lateral expansion (well before fNR) that is propagat- 
ing in a stellar wind external medium, which does poorly on 
point ( ii) above. 

We have found that for a jet with little lateral expansion be- 
fore fNR (our model 2), the jet would become roughly spher- 
ical only long after fNR (see Eq. |6land the discussion around 
it). This introduces a fast growth in the image size near fNR 
for on-axis observers with (?obs < 6*0 (see upper panel of Fig- 
ure |3} that detect the prompt gamma-ray emission (as in the 
case of GRB 030329). For an observer at ^obs ~ 90° as would 
typically be the case for GRBs that might be found in nearby 
SNe Ib/c, months to years after flie S N (Granot & Loeb 2003} 
'Paczvriski 2001: Ra mirez-Ruiz & M adau 2004), this causes a 
very slow increase in the image size near fNR (see lower panel 
o f Figure E). 

lOren. N^ar & PiranI j2004ft have considered a jet with no 
lateral spreading, even at f > fNR, and concluded that it can be 
ruled out for a uniform external density (k = 0) since it gives 
fl = 1/4 at fy < f < fNR which is inconsistent with observations 
[recall that in |6lwe have found that a = O.TIil^ j (I ct) between 
25 and 83 days]. In our analysis we have argued that physi- 
cally one expects lateral spreading to start around fNR, even if 
it is negligible at f <C fNR- We have shown that with this more 
realistic assumption for the jet dynamics (our model 2) the 
temporal evolution of the image size for a uniform external 
density (k = 0) is consistent with observations. 

The formalism developed in this paper would be 
useful for the analysis of future radio imaging of 
nearby GRB afterglows. The forthcoming Swift satellite 
(http://swift.gsfc.nasa.g0v/1 is likely to discover new GRBs 
at low redshifts. Follow-up imaging of their radio jets will 
constrain their physical properties and reveal whether the con- 
clusions we derived for GRB 030329 apply more generally to 
other relativistic explosions. 
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APPENDIX 

THE ANGLE 0* ON THE EQUAL ARRIVAL TIME SURFACE WHERE r± IS LOCATED 

The time at which a photon emitted at a lab frame time fiab and at spherical coordinates (r, 9, (j)) reaches the observer is given 



by 



r=flab-(^/c)C0s6' 



(Al) 



and shall be referred to as the observed time, where for convenience the direction to the observer was chosen to be along the 
z-axis (i.e. at 6* = 0). Let the location of a spherically symmetric shock front (or any other emitting surface for that matter) be 
described by r = ^(fiab) and that of an axially symmetric shock front by r = R(tiab,9). We shall now calculate the angle 9^ on the 
equal arrival time surface (which is defined by f = const) where R± = ma\(Rsm9) is located. At this point on the equal arrival 
time surface we have 



= 



/'dRsini. 
[ 89 



/dRsm9 
V 89 



8Rsm9 
8t\.. 



ab 



^\ =R(cos9 + Resm9) + f3rCsm9(^ 



where we use the notions (8R/8ti^\,)g = (3rC and Rg = {8\viR / 89)^.^^^. From Eq. ( lAU we have 

dt' 



= 



^ =^(sin0-^ecos0) + (l-/?,cos0)(^ 



(A2) 



(A3) 
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Table 1 . Comparing the values of Physical Parameters derived from dif- 
ferent Observables 



external 
density 


physical 
parameter 


observables i „ i„i i maior source uncertain by 

model 1 model 2 J -' 
being used of uncertamty a factor ot 


— n 
K — u 




F,(10d) 0.024u;%^ 0-401^,^^3* oc ~ 10-100 
Rx(0 0.8 5 oc r'', in model 1 (2) ~ 10 (~ 5) 


k = 2 


Esi/A, 


F,(10d) 03lu;'/^ 0.lOu;\{^ oci^a'''ym'''Ft'X^ -2-3 
i?x(0 2 1.2 oc s'^^''' in model 1 (2) ~ 5 (~ 2) 


k = 


r(25 d) 


„ 2.4 1.5 Rx(obs) & IX^',) ~1.3 
^ 2.1 2.4 Rx(obs)& jet model ~ 1.1 -1.2 


k = 2 


r(25 d) 


^ 2.4 2.4 «x(obs)&r(6»,) ~ 1.3 
^ 2.6 2.8 «x(obs)& jet model ~ 1.1 -1.2 


k = 




^ ^ 0.023 0.10 «x &r(6',)inEq.|8| ~ 10 
±. i'<i^.,» Q Q Q24 Rx & r(Fig.|6l in Eq.m -^5-10 

Fi,(10d) 0.241^^(3 0.0781/^1'*, model & value of ~3 


k = 2 




^ ^ 0.023 0.023 «x &r(e*)inEq.|8| ~ 10 
±. •^<y,„ Q Q20 0.017 Rx & r(Fig.|6j in Eq.lHl -.5-10 

F^(IOd) 0.361/^(3 0.171^^(3 model & value of ~3 



Note. — Estimates for the physical parameters of GRB 030329 derived from different observable quantities for different models of the jet lateral expansion. 
The value ofE/A is estimated from the spectrum at 10 days (upper line) and from the fit to the observed image size (lower line). The value of r(25 d) is evaluated 
both as T(9f ) according to ^(upper line) and from the fit to the observed image size (lower line). The value of in first two lines is evaluated first using Eg. 151 
with the values of r(25 d) from the con'esponding lines. In the third line the value of te is from the spectrum at 10 days (third line). 



SO that 

'9ftab\ R fRecosO-sinO 



Substituting Eq. (IA4> into Eq. iA2\ we obtain 



89 J , c \ l-/?,cos6i 



(A4) 



1 / dR\ sine /dR\ 

For a spherically symmetric shock Re = and cos 9^ = I3r{9*) = /?*, where in this case /3, is the shock velocity at the point on the 
equal arrival time surface were 9 = 9^, and R±^ is located. For a shock with axial symmetry we have 



P,-ReJ\-l3} + Rl 

cos 9, = L , (A6) 

\+Rl 



and 



I3, = I3,\\+R' (A7) 



where /3* is the shock velocity component normal to the shock fro nt in the rest frame of the upstream medium, which is the one 
that enters into the shock jump conditions jKumar & Granoi2003l) . 

THE APPARENT VELOCITY 

The apparent velocity, /3ap = {( \+z)/c\{dR±/dt), for a point source moving with a velocity /3 at an angle 9 from our line-of-sight 

is 

/3ap=/f^. (Bl) 

1-pcos^ 

For a spherical shock front movin g at a constant velocity /3sh, Ri_ is located at a constant angle 6'* which satisfies cos6'» = /?* = 
/3sh = const (according to Eq. IA5> so that the apparent velocity of the edge of the observed image is simply given by substituting 
COS0, = /3, in Eq. (IB1> . This gives 

/3,p = r,/3, = ^r2-i . (B2) 

We shall now show that this result holds for any spherically symmetric shock front. At t + dt we have 

9*(t + dt) = 9^(t) + d9^ , PJt + dt) = PJt) + dP:, , ^^/3, = <icos6', cx <if , (B3) 
and since Eq. (IB2> holds for a sphere moving at a constant velocity, we have 

[Rsm9](t + dt,9^) = R^(t) + T4t)P^(t)cdt + 0{dt^) . (B4) 
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Now, since is located where (dR sm9/d9)i = then 

R±(t + dt) = {R sin 6*] (t + dt,9^ + dO^) = {R sin OMt + dt ,9^) + 0{dt^) = R±(t) + r^(t)(34t)cdt + 0{dt^) , (B5) 



and therefore Eq. iB2\ holds for any spherically symmetric shock front. 
Finally, for an axially symmetric shock front, we obtain based on similar considerations as in the spherical case 



/3rSin6'» 
1-/3, cos 61* 



where 9^ and /3, are given by Eqs. IA6l and lX71 respectively. 

SOLVING FOR THE PHYSICAL PARAMETERS FROM A 'SNAPSHOT' SPECTRUM AT t > tj 

The most common ordering of the spectral break frequencies at tj <t < Jnr is < I'm < Vc, for which we obtain 



Vm=l35 X 10'^ (;'-0.67)(l +z)e'ey'«o'^%Tfd"ays Hz , (C2) 

Vc=l.l5 X 10'^(p-0.46)e-''^''(l +zr^fB'^nl^'^E';l''^{\ +Y)-^ Hz , (C3) 
i^..max = 1 3 1 (p + 0. 14)( 1 + zfefn'J^Elihll^fi-j^^^ mJy , (C4) 



for a uniform external medium {k = 0), and 



= 3.85 X l^ \^Xop'^2fl^ ^' ^zr'^'^>fAl^%'A2i' Hz , (C5) 



i.„,= 1.05 X 10'^ ( ^) (p-0.6m+z)e^/JW^E,iq^^^ Hz , (C6) 



p-1 



iy, = l.l5 X 10"(3.45-p)e°'*^''(l+zr'e/^A*^/^£5i(l+Fr2Hz , (C7) 

^2 1/2.1/2,7.. ,-1 r.-2 



/^,..max = 201(p + 0.12)(l+z)"e;,'^A;'^£5ifdays£>Z,28 mJy , (C8) 

for a stellar wind environment {k = 2), where F is the Compton y-parameter. A* = A/(5 x 10" grcm"'), fdays = t/(l day), eg 
is the fraction of the internal energy behind the shock in the magnetic field, and = 2/(10' x the c.g.s. units of Q). The 
emission depends only on the true energy in the jet, E, and does not depend on its initial half-opening angle 9q, since at f > tj 
(or equivalently when F drope below 9q^) the dynamics become independen t of i.e. the je t begins to expand sideways 
exponentially with radius in a self similar manner that is independent of 9o JGranot et al J 120021) . Solving the above sets of 
equations for the physical parameters yields 



3/2 

■/o(2.2) -'''^ VlmJy 



Mp) ,-5/6, -5/12,, 1/4 / ^t/.max \ .-1/2(1 +^-)-2n3 _n + y^l/2 

3/4 



£,so,52=o.io4^^z.;-v;-3-<- Q;:(i+zrDi,,ii+Yy'^ , (C9) 



17 Pi r,T2(^ SO.e(p) 5/12 5/24 3/8 / -^V.max \ ,5/4,, , N-lr,3/2,, , v\3/4 ^r-im 

£51 =0.0136 ^^^^^ ^^ J^Jg <13<14 (j-^ ) fdays(l+2) ^l;28(1+5^) > (^10) 

nmi/i Snip) 25/6 25/12 3/4 /-^V.max \ .7/2,, , s5r>-3 /1 , v\3/2 

«0 = 0.0714-^Z^„9' Z/„,(3 l/^,;i4 I ) fdays(l+2) ^L,28(l+5') ' , (Cll) 



1/2 

■gO,B(2.2)''«'y -m,13-cJ4 yi^jy 



gnMP) .-5/2.-5/4.-5/4 ^Wnax^ '-f/^2(l+z)-3DL28(l+Fr'/^ 

1/2 



., = 2.42^:^^.-' .-'3.;-; Q/;(i + ,)-z),,28(l + Fr/^ , (C12) 



n occ 80,e(p) 5/6 11/12 1/4 / ^i/.maA 1 .J/^ /1 , \ r->-i /1 . vNi// //-inN 

'^ = ^-^^^ia2:2)"«.^ "-'3 "-'4 [u^ ) W(l + ^)^"8(1+J') , (C13) 

for a uniform density, where fo(p) = e^^-^'^i'ip- iy/^(3p + 2)-^^^(p-0.67f/^\p-0A6T^/'^(p + 0.U)-^/^, goMp) = e°-^^^P(p- 
l)-^/^(3p + 2y/\p-0.67)-^/-\p-0A6)-^/\p + 0.U)-^^\ g„(p) = e^-^'^''(p-l)-^'H3p + 2y/\p-0.67)-^^'^Hp-0A6)-^/\p + 
OA4)y^,goAp) = e-^'^^P(p-iy'H3p+2)-y\p-0.6iy/\p-OA6y^\p+OA4)-^/\goAp) = e^^^^ 
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0.67) "/'^(;?-0.46) '/"^(p + O.U)'/^. For a stellar wind environment we find 

c n n^;^/l f^^P^ -5/6 -5/12 1/4 /■F^.max\ -1/2,, , s-2t-i3 /i , v\1/2 rr^i a\ 

£iso,52 = 0.0674-^r^J^J I'^^i'a i^.V ( ) ^£,28(1+^^)' , (C14) 

c AAi*;i 82.e(p) 5/12 5/24 3/8 /-PI, n,ax V'''^ .5/4 ,, , N-lr,3/2/i , t/x3/4 //^icn 
£51 =0.0161 ^^^^^ <13'^..14 ( Yml^ j fd.ys(l+^) ^L'28(l+i^) . (C15) 

^*=0.0262^^^.f./Vf/f3^;5fi,,(l+z)(l+y) , (C16) 

1 /2 

n (:on S2,b(p) -5/2 -5/4 -5/4 /fiy.max \ .-5/2,, , ._3„ , v\-5/2 ^i^-in\ 

^B = OMO—^iyJiy„, {^i^J^ i TmF^ ) W.s*^l+^^ 0^.28(1 + 7) ' , (C17) 

n co-r ^2,e(;') 5/6 11/12 I/4/fi/.max\ 3/2,, , /1 , v\l/2 o\ 

£, = 0.526 ^^^^^ ^^ lyjg lyj^ ly/^^ i TmF^ ) 'days(l +^)^l,28(1 ' , (C18) 

where /2(/?) = e-°-^^^P{p- iy/\3p + 2)-^/\p-0.69)^'^H^-'^5- pT^^\p + 0.12)-^/^, g2.,E(p) = e-°-^^'>P{p- l)-^^\3p + 2y^\p- 
Q.69)^^\3A5 - p)^^\p +0.12)^^ gAip) = e-^^-^^^P{p- l)-\3p + 2)(p-Q.69)-^/^(3A5- p)-^^^, gi.sip) = e°-^^^P(p - l)^^\3p + 
2r3/2(p-0.69)5/4(3.45-p)5/4(p+0.12ri/2, g2.e{p) = e-^^-^^^P(p-2)-Hp-iy/H3p+2y^\p-Q.69)-^^/^H3A5-p)-^/\p+0.12y/^. 

As was pointed out by Sari & Esin (2001), the expressions for the physical parameters that are derived from the instantaneous 
('snapshot') spectrum do not depend on the external density profile (i.e. on the value of k in our case), up to factors of order 
unity. This is because the instantaneous spectrum samples only the instantaneous external density just in front of the afterglow 
shock, next('")- The expression for the external density n for a uniform medium (k = 0) represents the density just in front of the 
shock for a general density profile that varies smoothly and gradually with radius, n < — > next('")^ where in our case «ext =Ar~^ /nip. 
However, for a non-uniform density next changes with radius and therefore with time. In our case, we assume the functional form 
of yisxtif) is known (i.e. we fix the value of k) and express the density normalization A as a function of the instantaneous values 
of the peak flux and break frequencies. 

We note that the expressions for the physical parameters at tj < t < Tnr are identical to those at f < tj. This is because we 
assume that the jet is uniform within a half-opening angle dj w F"', and therefore its emission is practically indistinguishable 
from that of a spherical blast wave with the same Lorentz factor F and radius R, or equivalently'" the same isotropic equivalent 
energy Eiso (which for a spherical blast wave is equal to the true energy, and for a model 1 jet is Eiso ~ (2 /dj)E « 2F^£') and 
observed time t (for the same values of next, and p). 

At / < tj, Eifio = const and is the more interesting physical quantity, while E in Eqs. (IC10> and (IC15> represents the energy within 
an angle of F"' around our line-of-sight which has no special physical significance at this stage. At tj < t < Tnr, however, the 
situation is reversed and E = const represents the true kinetic energy of the jet, and is therefore of great interest, while £150 ~ 2F^£' 
decreases with time and is no longer a very interesting physical quantity. 

For model 2, the jet continues to evolve as if it were part of a spherical blast wave with the same Eiso until fNR(£iso), and 
Eiso ~ (2/ Ol)E = const. Therefore, the emission at tj < t < Tnr is the same as from a spherical blast wave with the same Eiso, 
except for the peak flux /v.max which is suppressed by a factor of ^ (^oF)^ w (f/fy)"'-^"'^'/*^"*^ Hence, the above equations for the 
physical parameters may still be used in this case with the substitution /v.max — > ^V.max(f / tjf^~'^^l'-'^'^K In addition to this, in order 
to obtain the true energy in the jet, the expression for E (Eqs. IC10l and lC15> should be multiplied by (f / f^)"(3~'^)/(4-*)^ which is the 
fraction of the area within an angle of F"' around the line-of-sight which is occupied by the jet. 
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